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Abstract 

Accurate cervical intra-epithelial neoplasia (CIN) lesion grading is needed for effective patient management. We applied 
computer-assisted scanning and analytic approaches to immuno-stained CIN lesion sections to more accurately delineate 
disease states and decipher cell proliferation impacts from HPV and smoking within individual epithelial layers. A patient 
cohort undergoing cervical screening was identified (n = 196) and biopsies of varying disease grades and with intact 
basement membranes and epithelial layers were obtained (n = 261). Specimens were sectioned, stained (Mibl), and scanned 
using a high-resolution imaging system. We achieved semi-automated delineation of proliferation status and epithelial cell 
layers using Otsu segmentation, manual image review, Voronoi tessellation, and immuno-staining. Data were interrogated 
against l<nown status for HPV infection, smoking, and disease grade. We observed increased cell proliferation and decreased 
epithelial thickness with increased disease grade (when analyzing the epithelium at full thickness). Analysis within individual 
cell layers showed a >50% increase in cell proliferation for CIN2 vs. CINl lesions in higher epithelial layers (with minimal 
differences seen in basal/parabasal layers). Higher rates of proliferation for HPV-positive vs. -negative cases were seen in 
epithelial layers beyond the basal/parabasal layers in normal and CINl tissues. Comparing smokers vs. non-smokers, we 
observed increased cell proliferation in parabasal (low and high grade lesions) and basal layers (high grade only). In sum, we 
report CIN grade-specific differences in cell proliferation within individual epithelial layers. We also show HPV and smoking 
impacts on cell layer-specific proliferation. Our findings yield insight into CIN progression biology and demonstrate that 
rigorous, semi-automated imaging of histopathological specimens may be applied to improve disease grading accuracy. 
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Introduction 

Predicting outcomes for cervical intra-epithelial neoplasia (CIN) 
lesions remains a complex challenge. Some lesions progress to later 
disease stages while others do not, meaning some patients 
experience risks and costs of treatment unnecessarily. Further, 
HPV infection status for normal and early CIN tissues may be 
insufficient for stratifying progression risk. New tests are needed to 
accurately stratify patients presenting with CIN and to reduce the 
number of women treated unnecessarily for high grade squamous 
intraepithelial lesions (HGSILs). 

Multiple biomarkers have been tested to identify CIN lesions 
with a high risk of progression. P53, pl6, and Ki67/Mibl are 
amongst the best accepted for patient management [1-5]. It is 
known that proportions of proliferating cells increase with 
dysplastic stage. Recently, combined Mibl and pl6 analysis 
separated HGSILs based on progression risk [6]. Validation and 



acceptance of any biomarker requires an understanding of the 
molecular role that marker plays in disease. Others have analyzed 
the capacity of Mib 1 expression to identify high risk lesions and 
assist in diagnosis of HGSIL. Some groups have developed 
algorithms to quantify the distribution of proliferating cells and 
have demonstrated the power of these quantitative features over 
conventional, subjective assessments [3,5,7-13]. 

It is well-accepted that smoking is a cofactor for development of 
CIN [14—27]. Diverse hypotheses attempt to explain the effects of 
smoking, however, while smoking is recognized as a CIN co- 
factor, the exact nature of interactions between smoking, HPV 
infection, and dysplasia remains unclear [18]. 

Herein, we report our analysis on the effects of HPV infection 
and smoking on cell proliferation for normal and neoplastic 
cervical epitheha. We sought to apply a rigorous semi-automated 
approach to quantify these effects in CIN lesions. To achieve this, 
we analyzed individual epithelial layers in a well-annotated. 
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thoroughly reviewed patient cohort. Through this, we have gained 
insights into the impact of these factors on cell behavior for 
different disease stages. This work provides a rationale for wider 
evaluation of a combined approach involving clinical features (e.g. 
HPV, smoking status) and automated analysis of protein 
expression in epithelial layers as a biomarker for managing CIN. 

Materials and Methods 

Sample collection 

Samples were chosen amongst 1850 patients (3735 biopsies) 
aged S:18 that were collected during a multi-center study to 
evaluate Mib 1 and p 1 6 staining as a means of improving diagnosis 
of HGSIL [6]. Enrolled patients were those from a diagnostic 
population (i.e. had previously had an abnormal Pap test result). 
Seeking a distribution of lesion types and a cohort sufficiendy large 
to power meaningful statistical analyses, we chose 196 patients for 
whom 453 biopsies were available. As this work sought to analyze 
the fuU thickness of tissue cell layers, we further restricted our 
sample cohort to those biopsy specimens with intact basement 
membranes and epithelial layers at pathology review (n = 261). 
Samples were chosen from those acquired at the British Columbia 
Cancer Agency (BCCA), one of four institutions in the wider 
study. All samples meeting the above criteria (n = 26 1) - including 
multiple samples from the same patient - vmderwent Mibl 
staining and quantitative imaging (described below). 

Patient samples 

Non-pregnant women &18 years were enrolled from 1999- 
2006. The study protocol was approved by the University of 
British Columbia/BC Cancer Agency Research Ethics Board 
(protocol #C02-0476) and written consent was obtained from all 
enrolled patients. Patient demographic features (Table 1) were 
obtained for all cases, with ethnicity categories defined based on 
US OMB criteria, age defined as the subject's age at the time of 
the clinic visit, and smoking status defined based on current 
smoking habits. Regarding disease grade, biopsies were taken from 
one or two colposcopically abnormal areas and from two 
colposcopically-normal areas. Consensus diagnoses, based on 
review by expert pathologists, were used as the guiding criterion 
standard (see following section) [28]. 



Consensus diagnoses 

Histopathological specimens were reviewed at BCCA by both 
WHO and Bethesda criteria. Initial reviews were performed by 
one of the gynecological pathologists on clinical duty (the "pool"). 
There were seven pathologists in the pool, which also included 
study pathologists. A second blinded review was performed by one 
of two study pathologists. If the first and second readings agreed 
exactly in the WHO system, third reviews were not performed. 
Discrepancy of two grades based on WHO criteria mandated a 
third blinded reading by both study pathologists. Kappa scores for 
pathology readings have been reported [28]. In sum, for all 
diagnoses ranging from normal through CIN lesions to cancer, 
kappa-values ranged from 0.40 to 0.80; those in the HGSIL/ 
cancer range were in the high 0.70-0.85 kappa-value range. While 
"biopsy diagnoses" described disease state for a given biopsy 
specimen, "patient diagnoses" described the highest grade disease 
detected by biopsy review in a patient (since multiple biopsies were 
obtained from each patient and different disease stages were 
sometimes identified). 

HPV testing 

Endocervical samples were collected and tested for HPV using 
the Hybrid-Capture II system. This assay detects both low-risk 
(HPV-6, -11, -42, -43, -44) and high-risk (HPV-16, -18, -31, -33, - 
35, -39, -45, -51, -52, -56, -58, -59, -68) HPV types. All samples 
were processed using funds from pubhc grants. 

Mibl staining 

Paraffin sections (4 |im) were mounted on slides, dried 
overnight (37°C), deparaSinised in xylenes, and rehydrated in 
alcohol. Next, endogenous peroxidase activity was blocked by 3% 
H2O2 in phosphate-bufiFered saline (PBS). Sections were immersed 
in sodium citrate buffer and microwave heated at power settings of 
1000 W for 2 min, then 160 W for 15 min. Before immuno- 
staining, slides were soaked in PBS. Sections were incubated with 
biotinylated swine anti-rabbit antibody (DAKO) at 1:100 for 30 
min. Visualization of the complex was realized using diamino- 
benxidine/H202 for 10 min (room temperature). Two washes in 
PBS were performed before counterstaining with Mayer's 
Hematoxylin. Sections were dehydrated using graded ethanol 
and mounted with a mixture of distyrene, plasticiser, and xylenes. 



Table 1. Dennographic features for study patients. 





Patient Demographics 




n = 196 


Median Age (range) 




32 (18-66) 


Ethnicity 




White 


143 




Black 


2 




Hispanic 


3 




Native American 


4 




Asian 


30 




Other 


14 


Smol<ing status 




Current 


55 




Former 


46 




Never 


95 



doi:l 0.1 371 /journal.pone.Ol 07088.t001 
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Figure 1. Spatial analysis of iVIibI -stained cervical biopsies. We present: (A) an unprocessed, stained cervical biopsy cross section; (B) basal 
and superficial membranes that were first manually delineated by a technician through simple thresholding followed by automated segmentation of 
all nuclei, which generated candidate cell's nuclear centers of gravity; (C) Mibi -positive nuclei that were identified manually by the technician (green 
dots) and Voronoi diagrams that were generated based on those centers of gravity; (D) the basal layer (bottom of panel) containing all nuclei whose 
associated Voronoi polygons intersected with the basal membrane; (E) layer 1 (the parabasal layer); and (F) successive layers that were incrementally 
calculated. 

doi:1 0.1 371/journal.pone.01 07088.g001 



Table 2. Results for Mibl -stained biopsy specimens (including stratification based on histopathological review of an individual 
biopsy specimen ["Biopsy Diagnosis"]; worst disease grade observed amongst multiple biopsies from a single patient [i.e. "Worst 
Patient Diagnosis"]; and infection status for High Risk HPV types [i.e. the 13 high risk types defined within the Hybrid Capture II 
system]). 





Biopsy Being 
Investigated Diagnosis 


Worst Patient 
Diagnosis 


Negative for Hybrid Capture 
11 HPV high risk type 


Positive for Hybrid Capture 
11 HPV high risk type 


Total Number of 
biopsies studied 


Normal 


Normal 


15 


24 


39 


Normal 


LSIL 


28 


43 


71 


Normal 


HSIL 


2 


49 


51 


LSIL 


LSIL 


11 


16 


27 


LSIL 


HSIL 


0 


8 


8 


HSIL 


HSIL 


2 


63 


65 


Total 




58 


203 


261 



doi:l 0.1 371 /journal.pone.01 07088.t002 
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Table 3. Mean value of epithelial thickness (|xm), number of layers, and percentage of proliferating cells according to 
histopathological diagnosis (with standard deviations represented in parentheses. 





Pathological 
Diagnosis 


Total Number 


Patient Age 


Thickness (|im) 


Number of 
Layers 


% Proliferating Cells 


Normal 


39 


37.6 (11.0) 


420.7 (216.0) 


9.9 (3.0) 


15.1 (8.6) 


CINl 


35 


32.9 (8.9) 


331.4 (196.0) 


9.6 (3.5) 


28.7 (17.2) 


CIN2 


31 


30.3 (6.5) 


334.2 (200.0) 


1 1 .0 (4.4) 


34.3 (15.1) 


CIN3 


34 


29.6 (7.5) 


239.0 (125.0) 


12.1 (5.4) 


54.7 (19.7) 



doi:l 0.1 371 /journal.pone.Ol 07088.t003 



Interpretation of immuno-histochemistry staining 

Getafics [29] , an in-house high-resolution image analysis system 
using a monochrome Charge Coupled Device (CCD) camera 
mounted on a light microscope and attached to an image analysis 
software package, acquired image data. The illumination wave- 
length was 600 ±5 nm, effective pixel sampling space within the 
sample plane was 0.34 |J,m , and the effective pixel sampling area 
was 0.116 |J,m^. In each case, areas corresponding to diagnostic 
areas selected by study pathologists were carefijUy located by 
expert technicians. Images of that region were then captured. 
Regions of interest (ROI) were interactively and carefully 
delineated. Both basement membranes and superficial apical 
surfaces were separately defined. 

A simple, automated thresholding method followed by an Otsu 
segmentation algorithm was applied to extract centers of gravity of 
all cell nuclei [30]. Each segmented image was revisited and 
manual correction was performed by cyto-technicians to identify 
nuclei "forgotten" by automated algorithms and to remove false- 
positive objects. Mibl -positive nuclei were interactively marked, 
with Mibl status verified under microscope when necessary. 



Coordinates of the two membranes and spatial coordinates of the 
nuclei with negative/positive Mibl staining results were then 
saved and used for mathematical analyses. Tangentially cut 
epithelia were removed from analysis. 

Voronoi diagrams, layers-based analyses, and Voronoi 
layer definition 

Voronoi cell layers were defined thusly: all cells whose 
corresponding Voronoi polygons touched the basement mem- 
brane were assigned to layer 0 (basal layer). Layer 1 (parabasal 
layer) consisted of those cells which were not in layer 0 and had a 
neighbor in layer 0. Layer 2 consisted of those cells which were not 
in layer 1 or layer 0 and had a neighbor in layer 1 . Higher number 
layers were defined similarly. Generally, layers from layer 2 to the 
penultimate layer observed were interpreted as "intermediate" 
epithelial layers, while final layers were termed the "superficial" 
epithelial layer. 

The percentage of Mibl -positive cells in each layer was 
calculated as the number of positive cells divided by the number 
of cells belonging in the layer. Voronoi diagrams and layers-based 
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Figure 2. Epithelial thickness of normal epithelium in patients with no abnormalities, patients with LGSIL (i.e. CIN1), and patients 
with HGSIL (i.e. CIN2/3). No variability due to age was observed. The mid-point represents the median value; boxes represent the 25"^ percentiles; 
and whiskers represent the 95th percentiles. Factorial ANOVA testing was followed by a post hoc Fisher LSD test for group-by-group comparisons, 
yielding: p = 0.97 for a normal vs. LGSIL comparison; p = 0.048 for a normal vs. HGSIL comparison; and p = 0.020 for a LGSIL vs. HGSIL comparison. 
doi:1 0.1 371/journal.pone.01 07088.g002 
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Figure 3. Proportion of Mib1 -positive cells (i.e. proliferating cells) in the four diagnostic groups. The mid-point represents the median 
value; boxes represent the 25th percentiles; and whiskers represent the 95"^ percentiles. Once again, factorial ANOVA testing was followed by a posf 
hoc Fisher LSD test for group-by-group comparisons, yielding: p = 0.0006 for a normal vs. CIN1 comparison; p = 0.017 for a normal CIN1 vs. CIN2 
comparison; and p = 0.00001 for a CIN2 vs. CIN3 comparison. 
doi:1 0.1 371 /journal.pone.01 07088.g003 



analyses were undertaken as follows: given a set of points S in a 
plane, a Voronoi tessellation for S is the partition of the plane 
which associates a polygon V(p) with each point p of S in such a 
way that all locations in V(p) are closer to p than to any other point 
in S [31]. Applied to centers of nuclei in the ROI, the Voronoi 
tessellation allowed each nucleus to be associated with a Voronoi 
polygon that could be viewed as its area of influence. In this way, 
the Voronoi diagram gives a very natural mathematical represen- 
tation of the epithelium. With this diagram, it is straightforward to 
analyze the distribution of nuclei within a given region. We used 
structural features based on the Voronoi diagram and on Mibl 
positivity to describe the spatial interrelations of Mib 1 -positive 
cells in conjunction with the tissue architecture. Given a finite set 
of points, the Voronoi polygons of those points on the outside 
edges of the set were unbounded (i.e. one vertex of those polygons 
wiU be located at infinity). The polygons corresponding to these 
points were called marginal polygons. For this reason, the 
intersection of a marginal polygon with the ROI (consisting of 
the basement membrane, the external membrane, and the lines 
joining the ends of these two membranes) was found. Each 
marginal polygon was replaced with the polygon resulting from 
the intersection of the marginal polygon with the ROI boundary. 
Additional details related to Voronoi tessellations are provided in 
Fig. 1. 

Statistical Analysis 

AU analyses were done with STATISTICAIO. All p-values were 
two-sided, with a two-sided p-value<0.05 considered significant. 
For comparisons between groups, ANOVA testing was done. 
When testing showed a statistical difference across groups. Least 
Significant Differences testing or planned comparisons were used 
to compare groups two-by-two. 



Results 

Two-hundred sixty-one biopsies obtained from this diagnostic 
population (i.e. cohort with a previous abnormal Pap test result) 
were stained with Mibl (Table 2). Based on worst observed patient 
diagnosis in each case, 62% (24/39) of all "normal" samples were 
HPV-positive. Also based on the worst observed patient diagnosis, 
61% (59/97) of Low-Grade squamous intraepithelial lesions 
(LGSILs, i.e. CINl lesions) were positive for oncogenic HPV, 
whereas 97% (122/126) of HGSILs (CIN2/3 lesions) were HPV- 
positive. Normal epithelia from cervices in which a HGSIL was 
found were on average thinner than epithelia from women in 
which no abnormalities were identified - and also from epithelia 
from women with LGSIL (data not shown). These differences were 
statistically significant (p = 0.04 and p = 0.02, respectively). There 
were no differences between the percentage of proliferating cells 
and the number of layers in these groups. Nevertheless, because of 
these differences in epithelium thickness, we restricted analyses to 
normal specimens from patients without abnormalities when 
describing normal tissues (Table 3). Hence, downstream analyses 
spanned i) all HGSIL biopsies, ii) all LGSIL biopsies, and iii) all 
"normal" biopsies from cervices of patients where no abnormality 
was detected in any biopsy. 

The mean and standard deviation of epithelial thickness, 
number of layers, and percentage of proliferating cells for four 
pathological grades were determined (Table 3). There was a 
consistent decrease in epithelial thickness from normal epithelium 
to CIN3 lesions. Epithelia of CINl lesions were on average 22% 
(89.3 |J.m) thinner than those of normal tissues. For CIN3 lesions, 
epithelia were 28% (95.2 |J.m) thinner than in CIN2 lesions. There 
was a statistically significant difference between normal and CIN3 
tissues (p<0.0001). Additional epithefial thickness comparisons 
based on disease grade were undertaken (Fig. 2). 
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Voronoi layers were calculated as described above and in Fig. 1 , 
The number of layers increased from normal to CIN3 tissues, with 
a statistically significant difference observed between CIN3 vs. 
normal specimens (p = 0.03) and between CIN3 vs. CINl lesions 
(p = 0.021) (Table 3). Fig. 3 compares the percentage of prolifer- 
ating cells vs. disease grade. A clear trend shows increased 
percentages of proliferating cells vs. increased CIN grade. These 
increases were statistically significant for normal vs. CINl tissues 
(p<0.001) and for CIN2 vs. CIN3 lesions (p<0.00001). 
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Layer-based cell proliferation and dysplasia 

In normal epithelia, >60% of cell proliferation occurred in the 
parabasal layer and layer 2 (Table 4). Proliferation in the upper 
layers dropped rapidly to 14%, 8%, and 4% in layers 3, 4, and 5, 
respectively - and was almost negligible for higher layers. 
Interestingly, we observed a near-linear increase in the percentage 
of proliferating cells in both basal and parabasal layers from 
normal epithelium to CINl, CINl to CIN2, and CIN2 to CIN3. 
The differentiation between CINl and CIN2 occurred in layer 3, 
where proliferating cells represented 50% of cells in CIN2 lesions 
vs. only 33% in CINl lesions. In higher layers, the proliferation 
profile of CINl lesions resembled the proliferation pattern of 
normal epithelial tissues. ContrarUy, CIN2 lesions exliibited much 
higher proliferation rates, with 30-40% more proliferating cells in 
each higher layer. 

We next returned to assessing proliferation on a layer-by-layer 
basis for the four diagnostic groups (Fig. 4A). The percentage of 
proliferating cells in the basal layer was much higher in CIN3 
lesions vs. normal epithelial tissues (30% vs. 6%) and almost 
identical between CINl and CIN2 lesions (—22%). A similar 
pattern was observed in the parabasal layer (layer 1): 60% of CIN3 
nuclei were Mibl-positive vs. only 33% for normal cases. CINl 
and CIN2 lesions exhibited similar percentages of Mib 1 -positive 
nuclei (~44%). For each of the three CIN groups, we noted that 
the percentage of Mib 1 -positive nuclei in the parabasal layer was 
approximately double the percentage in the basal layer: 30% vs. 
60% for CIN3 cases and 20% vs. 40% for CINl and CIN2 cases. 
In normal epithelia, the difference was larger: 33% of Mibl- 
positive nuclei were in the parabasal layer vs. 7% in the basal 
layer. 

This relative proliferation pattern changed in upper layers. The 
percentage of Mib 1 -positive nuclei for CINl cases decreased by 
~10% from layer 1 to 2 but increased by ~6% from layer 1 to 2 
in CIN2 lesions. The decrease in the percentage of Mib 1 -positive 
nuclei in CIN2 and CIN 1 lesions was almost identical from layer 2 
to layer 6, with an 18% decrease in CIN2 lesions and 14% 
decrease in CINl lesions, respectively. 
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Layer-based cell proliferation and dysplasia and HPV 
infection 

We noted several HPV impacts on proliferation patterns 
(Fig. 4B). This analysis was restricted to normal tissues and 
CINl lesions (too few HPV-negative HGSILs were available for 
meaningful analyses related to infection). We observed similar 
stratifications of proliferation profiles for normal and CIN 1 groups 
based on HPV status; for both normal and CINl cases, a higher 
percentage of Mib 1 -positive nuclei were noted in layer 3 for HPV- 
positive cases vs. HPV-negative ones. The layer proliferation 
profile of CINl HPV-positive tissues resembled the layer 
proliferation profile of the CIN2 lesions. (Proliferation profiles 
for women £40 vs. >40 years old showed near identical results, 
confirming absence of an age effect on cell proliferation rates [data 
not shown]). 
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Figure 4. Average percentage of proliferating ceils in the first layers of cervical epithelium. In all panels, layers 0 and 1 represent the 
basal and parabasal layers, respectively. A) Percentage of proliferating cells in the first seven layers of cervical epithelium as assessed for each 
pathology grade. B) Percentage of proliferating cells in the first seven layers of cervical epithelium for normal and CIN1 tissues based on high-risk HPV 
type infection status. C) Percentage of proliferating cells in the first six layers of epithelium of non-smokers, where LGSILs represent CIN1 lesions and 
HGSILs represent CIN2/3 lesions. D) Same representation as C), only results presented are for patients with smoking histories. 
doi:1 0.1 371/journal.pone.01 07088.g004 



Layer-based cell Proliferation and dysplasia and smoking 
status 

Fig. 5 shows the impact of smoking status on cell proliferation 
rates when the epithelium was analyzed at full thickness. Limited 
separation was observed between non-diseased and LGSIL tissues. 
The difference for HGSILs between smokers vs. non-smokers was 
not statistically significant. Fig. 4C and D shows proliferation 
profiles for the first 6 layers of normal, LGSIL, and HGSIL cases 
for non-smokers and smokers, respectively. Overall cell prolifer- 
ation was detected as higher for smokers for all layers. The most 
marked difference between smokers vs. non-smokers occurred in 
the basal layer of HGSILs, where the percentage of proliferating 
cells doubled for smokers (19% vs. 40%, p<0.001). 

Discussion 

Novel methods described herein allowed quantification of cell 
proliferation on a per-epithelial-layer basis (Fig. 1). This yielded 
insight into the spatial effects of smoking, HPV infection, and 
dysplastic processes in cell proliferation for CIN lesions. Combined 
analysis of all epithelial layers revealed an increase in cell 
proliferation with progression through dysplastic stages, agreeing 
with previous work (Fig. 3, Table 4) [1,12]. On average, the 



percentage of proliferating cells increased with progression: from 
16% in normal epithelia to 27% in CINl, 36% in CIN2, and 55% 
in CIN3 (Fig. 3). This increase in proliferating cells was 
accompanied by a decrease in epithelial thickness through loss of 
differentiated cells that were replaced by newly dividing, smaller 
cells (Fig. 2, Table 3). Consequently, cell density increased 
together with the number of cell layers. While proliferation rates 
for various CIN stages were comparable to those of previous 
reports, we observed an elevated proliferation rate in normal 
tissues [32] . This result could have been influenced by whether or 
not "histologically normal tissues" were taken from both diseased 
and non-diseased patients and also by general variability in disease 
state classification for LGSILs vs. normal tissues. 

One key observation was that differences in proliferation rates 
between CINl and CIN2 lesions occurred in specific epithelial 
layers (Fig. 4A). While proliferation rates for these lesions were 
nearly identical for basal and parabasal layers, proliferation rates 
in CIN2 lesions were £50% higher in all higher layers. Whereas 
combined analysis of all epithelial layers of CIN 1 and CIN2 lesions 
is not sufficient to discriminate these disease stages, this result 
suggests that honing analysis within epithelial layers above basal/ 
parabasal layers may yield more accurate clinical grading (which 
would impact estimation of disease progression risk). Increased 
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Figure 5. Percentage of proliferating cells in the full thickness of the epithelium for normal epithelium, LGSILs, and HGSILs for non- 
smokers and smokers. Factorial ANOV A/posf hoc Fisher LSD testing revealed the following results when comparing results from non-smokers vs. 
smokers: for cases negative for disease, p = 0.75; for LGSILs, p = 0.94; and for HGSILs, p = 0.049. 
doi:1 0.1 371 /journal.pone.01 07088.g005 



proliferation rates in higher epitlielial layers for CIN2 lesions (vs. 
CINl lesions) may be driven by the impacts of greater viral 
genome integration (e.g. disruption of viral E2 expression), which 
is associated with disease progression [33]. Interestingly, while no 
difference in basal layer proliferation rates was observed for CINl 
vs. CIN2 lesions, we noted a marked increase in basal layer 
proliferation rates for CIN3 lesions. This increase may be 
attributable to increased disease signalling in higher epithelial 
layers, a feedback mechanism that may further disease progres- 
sion. 

We also observed HPV-driven proliferation rate differences 
within a given disease stage. Moreover, these differences were 
detected in specific epithelial layers. Given high prevalence of 
HPV infection in advanced disease stages, we were only able to 
derive statistically meaningful results related to HPV status for 
normal and CINl tissues. For both these tissue types, we noted 
higher rates of proliferation for HPV-positive cases relative to 
HPV-negative cases in epithelial layers above basal and parabasal 
layers (Fig. 4B). These data suggest that analysis of individual 
epithelial layers may be an effective, independent means of 
determining infection status/ disease stage - and that assessing cell 
proliferation rates in each epithelial layer may be a useful tool for 
guiding disease management decisions. These results also suggest 
that CINl HPV-negative lesions may be driven by a hyperplastic 
response that up-regulates proliferation of basal cells - and that 
daughter cells subsequendy respond in the same fashion (which 
would suggest that they in turn foUow the same control sequence). 
In HPV-positive cases, it may be that HPV infection disrupts this 
control sequence, resulting in a different pattern of proliferation 
(vs. non-infected tissues). The molecular mechanisms by which 
higher proliferation rates in HPV-positive normal and CINl 
lesions occur only in epithelial levels beyond the basal and 
parabasal layers remains unknown. Follow-up analysis of a larger 
panel of HPV-positive/-negative CIN2 lesions is needed to 



determine whether proliferation in higher epithelial layers can 
delineate more advanced disease stages. 

Finally, we noticed marked smoking-mediated differences in the 
distribution of proliferating cells in CIN lesions. When epithelia 
were analyzed at fuU thickness, negligible separation in the 
percentage of proliferating cells was observed for normal and 
LGSIL tissues from smoking vs. non-smoking patients (a small 
degree of separation was observed for HGSILs, with smoking cases 
exhibiting somewhat higher percentages of proliferating cells) 
(Fig. 5). Nonetheless, separation into individual epithelial layers 
revealed doubled proliferation rates in the basal layer of HGSILs 
from smokers vs. non-smokers (along with increased proliferation 
rates in the parabasal layer) (Fig. 4C, D). This may be indicative of 
smoking-related factors driving cells out of senescence, which 
could in turn result in a larger pool of cells harboring molecular 
alterations, increasing the likelihood that a given cell in the tissue 
field win mutate further and progress to later stage disease. 
Interestingly, we also noted significantiy higher proliferation rates 
in the parabasal layer of smokers vs. non-smokers for LGSILs but 
no significant difference in basal layer proliferation for these cases. 
The striking difference in basal layer proliferation rates for 
HGSILs may be attributable to inflammation effects due to 
smoking, which can in general be associated with increased cell 
proliferation in other tissue types [34-36]. The absence of this 
smoking-driven phenomenon in LGSILs may suggest that the 
degree of inflammation in these cases is insufficient to trigger 
higher proliferation rates. Our finding of no major differences in 
normal tissue cell proliferation rates for smokers vs. non-smokers 
agrees with previous reports [37]. Taken together, our data suggest 
a strong effect of smoking on cell proliferation that manifests 
independently of disease stage and most prominently in basal and 
parabasal epithelial layers. 



PLOS ONE I www.plosone.org 



8 



September 2014 | Volume 9 | Issue 9 | el 07088 



Variations in Cell Proliferation within Cervical Epithelial Layers 



Conclusion 

Our analyses of cell proliferation in normal cervical and CIN 
tissues have identified epithelial layer-specific features that 
delineate disease grades. Our results also show that proliferation 
within cell layers varies depending on HPV infection and smoking 
status. We have demonstrated the utility of computer-assisted 
histopathological slide analysis for identifying novel cell features 
associated with disease grade. More accurate clinical diagnoses 
based on this technique could improve patient management by 
facilitating more accurate estimations of disease progression risk. 
This work bolsters the argument that research into HPV and 
smoking impacts on cell proliferation must account for spatial 
differences in proliferation based on individual epithelial layers 
(instead of relying on analyses of the entire epithelium at its fuU 
thickness). Independent analyses using the approaches described 
here are needed to determine whether cell proliferation in 
epithelial layers beyond the basal and parabasal layers has utility 
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